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Femtosecond Photodissociation of Molecules Facilitated by Noise 
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We investigate the dynamics of diatomic molecules subjected to both a femtosecond mid-infrared 
laser pulse and Gaussian white noise. The stochastic Schrodinger equation with a Morse potential is 
used to describe the molecular vibrations under noise and the laser pulse. For weak laser intensity, 
well below the dissociation threshold, it is shown that one can find an optimum amount of noise 
that leads to a dramatic enhancement of the dissociation probability. The enhancement landscape 
which is shown as a function of both the noise and the laser strength, exhibits a global maximum. 
A frequency-resolved gain profile is recorded with a pump-probe set-up which is experimentally 
realizable. With this profile we identify the linear and nonlinear multiphoton processes created 
by the interplay between laser and noise and assess their relative contribution to the dissociation 
enhancement. 

PACS numbers: 02.50.Ey, 33.80.Gj, 42.50.Hz 



I. INTRODUCTION 

The efficient dissociation of a molecular bond by the 
application of ultrashort laser pulses is a long sought after 
goal in photochemistry [l], [2j. However, to achieve sig- 
nificant molecular dissociation with coherent ultrashort 
laser pulses is difficult due to the anharmonicity of the 
molecular vibrations 0] . Very high intensities of the laser 
field are therefore required to produce significant dissoci- 
ation, despite the lower value of the dissociation energy 
compared to the ionization potential. For such strong 
laser fields molecular tunnel ionization often masks the 
dissociation process [1, @ . 

Alternative strategies have been proposed to enhance 
the dissociation yield for weaker laser pulses. For in- 
stance, the approach using the linearly (or circularly) 
polarized chirped laser pulses have been proposed. Here, 
the frequency is designed to match the vibrational lad- 
der of a specific molecule, thereby facilitating the lad- 
der climbing of the molecule and eventually its dissocia- 
tion These ideas have been confirmed in the exper- 
iments, using chirped femtosecond mid-infrared (MIR) 
laser pulses. A large population transfer to higher vibra- 
tional levels have been achieved with diatomic molecules 
(HCl, NO etc) as well as in some polyatomic molecules 
(Cr(CO) 6 ) 0,lEH- Furthermore, the optimally tailored 
pulses by a closed loop control have also been used quite 
recently to dissociate a specific molecular bond in a poly- 
atomic molecule; it was possible to control in particular 
the branching ratio between two possible reaction chan- 
nels [Io|. 

Following a completely different approach to achieve 
the molecular dissociation, it has been proposed recently 
that one can instead use an incoherent random kicks 
[Tl| . In the same spirit, there has been a growing inter- 
est in understanding the role of random fluctuations on 
the quantum dynamics of atomic and molecular systems 
0, EL Ell- However, the molecular dissociation pro- 
cesses due to the simultaneous application of a laser and 



a noise source has so far not yet been considered. In do- 
ing so, we expect the coherent laser pulse to become more 
efficient when submerged in a small amount of noise. It 
is indeed well known that in some classical as well as 
quantum nonlinear systems (e.g. bistable systems), the 
noise can lower the nonlinear threshold via the stochastic 
resonance (SR) effect [HEl. 

It is also worth mention- 
ing here that SR also exists in some nonlinear chemical 
reactions [HI, [l6[ , but has not been demonstrated in the 
crucial process of molecular (photo)dissociation on a fem- 
tosecond time-scale yet. 

Here, we propose an approach to control the molecular 
dissociation by irradiating diatomic molecules by a fem- 
tosecond MIR laser pulse in the presence of noise. We 
show how and to what extent the presence of noise can 
reduce the dissociation threshold of diatomic molecules 
due to its nonlinear interaction with the laser pulse. This 
approach can be applicable to both polar and nonpolar 
molecules and is readily possible to use in experiments. 

We begin the main text of this work in Sec. II with 
a brief description of the molecular Hamiltonian includ- 
ing a coherent laser pulse and a stochastic term. In 
Sec. Ill we describe the fast-Fourier transform (FFT) 
split-operator method to solve the stochastic Schrodinger 
equation with absorbing boundary conditions, and focus 
on the relevant physical observable of interest. Our re- 
sults for the role of noise in the molecular photodisso- 
ciation, and other interesting features of the effect are 
discussed in Sec. IV. Finally, Sec. V concludes the paper 
with some perspectives. 



II. THE THEORETICAL APPROACH 

The total molecular Hamiltonian for the diatomic 
molecule under time-dependent external fields is given 
by (atomic units are used unless stated otherwise) 



H(t) =H - atf(t) - xF(t), 



(1) 
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FIG. 1: The dissociation probability Pl as function of the 
laser peak amplitude Fo over long time t (t » T p ). Here, 
w = 0.007 a.u., 8 = 0.0, and T p = Wtt/lj. The arrow indi- 
cates the value of Fo = 0.04 a.u. for which Pl ~ is far 
less than the dissociation threshold. The inset shows the cor- 
responding 15 optical cycle laser pulse F(t) (light solid line) 
and the schematic pictures of a diatomic molecule HF (two 
nuclei separated by R, with Ro, being the equilibrium inter- 
nuclear distance) and the Morse potential V(x) (thick solid 
line), with x = R — Ro- 



where Hq = 2— + V(x) describes vibrational motion of 
the molecule with reduced mass m in the Morse potential 

m 



V(x) = -D e +D e (l- exp(-px)y , 



(2) 



with well depth D e and length scale 0. The eigenenergies 
E n of the Morse oscillator Hq are given by 

E„ = Hu) e (n+l/2)[l-B(n+l/2)/2], < n < [j] (3) 

where uj e is the harmonic frequency, n& = [j] + 1 is the 
number of bound states with 

j = l/B-1/2, B = hf3(2mD e )- 1/2 , hw e = 2BD e . (4) 

The laser field is a nonresonant MIR femtosecond 
pulse, 



F(t) = f{t)F sm(ujt + 5) , 



(5) 



where u> denotes the angular frequency, and the effective 
peak amplitude Fo as well as the noise amplitude £(t ) 
(see below) contain the molecular dipole gradient pT| . 
We choose a smooth pulse envelope f(t) of the form 

f(t) = sin 2 (Trf/Tp), (6) 

where T p is the pulse duration of typically 15 optical 
cycles (see inset of Fig. 1). 

The noise term £(i) is a zero-mean Gaussian white 
noise having the following properties, 



(£(<)) = o, 

where D is the noise intensity [l9[ . 



(7) 
(8) 



III. TIME EVOLUTION AND OBSERVABLES 

A. Quantum stochastic propagation 

The time-evolution under Eq. |T]) is non-deterministic 
due to the stochastic nature of the Hamiltonian. A so- 
lution taking into account the external fields perturba- 
tively is not possible, because the fields of interest are 
so strong that they influence the dynamics substantially. 
Hence, we have to solve the full time-dependent stochas- 
tic Schrodinger equation, 



(9) 



The solution is accomplished by averaging over a suffi- 
cient number of deterministic solutions under different 
specific realizations r of the noise. 

For a given realization r, the solution of the stochas- 
tic Schrodinger equation amounts to solve the standard 
time-dependent Schrodinger equation, 



IVVW) =U r {t,t )\ll)(to)), 



(10) 



starting from an initial state |^(fo)) at time io- The 
stochastic evolution operator U r (t,t n ) is written as a 
product of operators propagating over a small time in- 
terval At, 



Ur(At) = exp 



-if{ +At H(x,t)dt 



(11) 



This short-time propagator can be formulated explicitely 
by evaluating the integral in the exponent in the 
Stratonovitch sense [lfj], leading to 



U r (At) = exp (-i(H (x) - xF{t))At) exp(ixV2DAt£ t ), 

(12) 

where £t is a Gaussian distributed random number of unit 
variance. The representation of Eq. (fT5|) illustrates how 
the coherent evolution of the system (with laser) aquires 
an additional random phase in the form of a momentum 
kick whose strength is rand om due to £t but scales with 
the amplitude of the noise over the time interval 

At. The short-time propagator given by Eq. (fT2")) can be 
easily realized numerically using the FFT split-operator 
approach with absorbing boundary conditions [201 ] . 



B. Dissociation probability 

We use the dissociation probability as a measure for 
the non-linear coupling of energy from the external fields 
to the molecule. It can be most conveniently expressed 
as a complement of the probability to find the molecule 
at time t in any of its bound states, which reads for a 
given realization r 



N b -1 



(13) 
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FIG. 2: The average dissociation probability of the molecule 
for noise only Pn (square) and for the simultaneous applica- 
tion of the laser pulse and the noise Pl+n (circle). The total 
number of realization N r = 100, and the laser peak amplitude 
is kept fixed at Fo = 0.04 a.u.. 



where for the initial state \if>o(t = 0)) we take the ground 
state of the Morse potential with energy F (Eq. 3). 
Typically we will average over a large number of realiza- 
tion (N r about 100) to calculate the average dissociation 
probability over long time too » T p , 
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FIG. 3: The factor of enhancement r\ versus the noise ampli- 
tude yfD for the 15 cycle long MIR pulse with Fo = 0.04 a.u. 
and uj = 0.007 a.u.. 



laser pulse unchanged, we add a small amount of noise 
to it. Fig. 2 shows the average dissociation probability 
Pl+n versus the noise amplitude \[D (squares). One 
can clearly see that by adding a small amount of noise, 
Pl+n rises rapidly, well beyond the increase induced by 
the noise alone (Pn, circles) which is applied for the same 
time T p as the laser pulse. 

The enhancement by noise can be quantified with the 
parameter 



In the following, we shall consider different physical sit- 
uations such as when the molecule is subjected to (i) the 
laser pulse alone (denoted by subsrcript L), (ii) the noise 
alone (N), or (iii) a combination of both fields (L+N). 



IV. RESULTS AND DISCUSSIONS 



A. Femtosecond photodissociation under noise 



P, 



L+N 



-Pa 



Pn 



(15) 



with Pq = Pl + Pn which is plotted in Fig. 3. By con- 
struction, 77 — ^ if one of the two fields, laser or noise, is 
much larger than the other. However, for a specfic noise 
amplitude, here at \pD = 0.02, a dramatic enhancement 
of the order of 10 5 is observed. 



Photodissociation of diatomic molecules by a femtosec- 
ond laser pulse is a highly non-linear processes, as can be 
seen in Fig. 1 for a 15 cycle MIR pulse at 10 = 0.007 a.u.. 
The dissociation probability Pj, versus the peak field am- 
plitude Fo exhibits a prototypical S-shape curve with a 
threshold like behavior. High intensities (i.e., Fq val- 
ues) are necessary for dissociation by a laser pulse alone. 
However, the addition of a weak amplitude noise to a 
subthreshold laser pulse may lead to substantial dissoci- 
ation. As we have shown before, such a combination is 
quite efficient for an otpimum choice of the noise ampli- 
tude in achieving ionization of an atom [Li, resembling 
the phenomenon of stochastic resonance [lfj . [ItJ • 

We choose a laser pulse of peak amplitude F = 
0.04 a.u., for which the dissociation probability is ex- 
tremely small, Pl < 10 -10 (see Fig. 1). Keeping this 



B. The net enhancement landscape 

Having kept the laser amplitude fixed so far to a sub- 
threshold value (Fo = 0-04 a.u.), we now investigate the 
enhancement factor as a function of both, laser and noise 
amplitudes. As can be seen in Fig. 4, ^(Vl?, Fq) shows a 
global maximum. Additional calculations have revealed 
that the shape of r)(>/D,Fo) is quite robust for differ- 
ent range of parameters, such as carrier-envelop phase 5, 
pulse duration T p , and even for different values of the 
laser frequency u> in the infrared (IR) regime. Hence, 
the enhancement landscape from Fig. 4 can be used as 
a guide to identify the range of enhancement possible in 
experiments on molecular photodissociation under noise. 
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FIG. 4: The enhancement landscape r/(V D, Fo) versus the 
peak pulse amplitude Fo and the noise amplitude \fD. Other 
parameters are kept fixed as in figure 3. Contour lines of equal 
enhancement are also shown on the (\/I?,Fo) plane. 



C. Frequency-resolved enhancement profile for HF 

In order to shed light on the mechanism underlying the 
stochastic enhancement of dissociation we take a closer 
look at the frequency resolved enhancement profile for 
the case of HF. This will reveal which frequencies pro- 
vide the gain within the braod band noise and will in turn 
clarify the non-linear photo processes responsible for the 
enhancement. 

For this purpose, we record the frequency-resolved 
molecular gain G{tu) (FRMG) employing a pump-probe 
type of setting without noise. We consider the HF 
molecule to be driven by the MIR laser pulse of Fo = 
0.04 a.u. and lo = 0.007 a.u. (pump pulse), and re- 
place the noise by a probe pulse of amplitude F p and of 
tunable frequency io v . The tunable probe pulse has the 
same duration and envelope as the fixed pump pulse. The 
probe amplitude (F p /F ~ 0.05) is such that it can only 
produce single photon transitions. For a HF molecule 
in its ground state, G(ui) is obtained by measuring the 
net energy absorbed as a function of oj p and is shown in 
Fig. 5(a). It exhibits several clearly resolved peaks which 
play a dominant role in the dissociation enhancement un- 
der broad band light. 

We have identified these peaks as originating from sin- 
gle or multiple photon processes for the level-structure 
of HF molecule. For convenience, these processes can 
be classified as follows. All photo processes that involve 
the lowest two levels of the HF are termed A-series, 
while the ones involving the next higher level (three lev- 
els) are termed B-series, and so on. Fig. 5(b) shows 
these excitation pathways for the given pump frequency 
u> = 0.007 a.u.. These basically involve n-photons of 
the pump (n = 0, ±1, ±2, ..) plus a single photon of the 
probe. Note that in addition to the resonant single pho- 
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FIG. 5: (a) The frequency-resolved gain G(ui) for the driven 
HF molecule measured using a simultaneous pump-probe set- 
ting without noise. G(u) for the bare molecule is also shown 
(black curve) for comparison, (b) Various single-photon and 
mutiphoton transitions corresponding to the peaks in G{oj) 
are classified as A-series, B-series, C-series that involve first 
two, three, four vibrational levels of HF, respectively. Thick 
red arrow denote the pump photon (uj = 0.007 a.u.) and thin 
black arrow is for tunable weak probe pulse. Pulse band- 
width is here negligible compared to the photon energies, (c) 
Enhancement curve for the chaotic light with perforated spec- 
trum of hole-widths 0.004 a.u. is compared with the broad- 
band light of BW=1.0 a.u.. 



ton transitions, second- and third-order processes play 
an important role. Thus, the diatomic molecule offers an 
interesting example where the role of higher order photo- 
processes can be clearly isolated, which is nontrivial in 
the case of an analogous atomic system (l3| . 

To highlight the contribution of the multiphoton tran- 
sitions in the enhancement, we have eliminated the first 
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four single-photon resonant frequencies from the noise 
spectrum. Such a filtered noise can be easily simulated 
by means of the chaotic light [Hj]. For a chosen hole- 
width of 0.008 a.u. centered at the resonant frequen- 
cies of HF, a drop of more than one order of magnitude 
in the maximum r\ is observed as shown in Fig. 5(c). 
Note that one is still left with a significant enhancement 
(r\ ~ 10 4 ), which can be clearly attributed to the higher 
order processes other than single photon resonant tran- 
sitions. Although we have considered specifically HF, 
the approach used here can be applied to other diatomic 
molecules. We have found similar enhancement features 
for HCl and H2 molecules within the framework of the 
Morse potential subjected to different MIR pulses. 

V. SUMMARY AND CONCLUSION 

We have shown that it is possible to significantly en- 
hance the dissociation probability of a diatomic molecule 
under a weak femtosecond MIR pulse by adding a small 
amount of noise to it. The net enhancement for a given 
laser pulse exhibits a single maximum as the noise level is 
varied. This maximum suggests an optimum noise ampli- 
tude to achieve the maximum dissociation, in analogy to 
the stochastic resonance phenomenon. The enhancement 



landscape, that is the enhancement as a function of laser 
amplitude and noise amplitude, exhibits a global maxi- 
mum. Analyzing the frequency resolved gain profile of 
the molecule using a pump-probe setting we have iden- 
tified the dominant frequencies and the corresponding 
physical processes activated by the noise. We conclude 
that in addition to single-photon resonant transitions, 
multiphoton transitions play a significant role. 

Similar effects are expected for other diatomic 
molecules if exposed to a combination of noise and laser 
pulse, where pulse duration and wavelength in the MIR 
regime may be chosen from a wide parameter range. Of 
particular interest in the future will be if the proposed 
nonlinear interplay can achieve selective bond breaking 
by enhancing dissociation of a particular bond in a poly- 
atomic molecule. 
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